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ABSTRACT 

We present a new analysis of an archived Chandra HETGS X-ray spectrum of the 
WR+O colliding wind binary 7 2 Velorum. The spectrum is dominated by emission 
lines from astrophysically abundant elements: Ne, Mg, Si, S and Fe. From a combina- 
tion of broad-band spectral analysis and an analysis of line flux ratios we infer a wide 
range of temperatures in the X-ray emitting plasma (^4-40 MK). As in the previously 
published analysis, we find the X-ray emission lines are essentially unshifted, with a 
mean FWHM of 1240 ± 30 km s -1 . Calculations of line profiles based on hydrody- 
namical simulations of the wind-wind collision predict lines that are blueshifted by 
a few hundred km s _1 . The lack of any observed shift in the lines may be evidence 
of a large shock-cone opening half-angle (> 85°), and we suggest this may be evi- 
dence of sudden radiative braking. From the R and G ratios measured from He-like 
forbidden- intercombination-resonance triplets we find evidence that the Mgxi emis- 
sion originates from hotter gas closer to the O star than the Sixin emission, which 
suggests that non-equilibrium ionization may be present. 
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1 INTRODUCTION 

7 2 Velorum (WR 11, HD 68273) is the closest known Wolf- 
Rayet (W R) star, at a _ffijPjparcos-determined distance of 
258^31 pc iSchaerer. Schmutz fe Grenorjll997^ . As such, it 
is a key system for increasing our understanding of X-ray 
emission from massive early-type stars. r y 2 Vel is a double- 
lined spectroscopic binary of spectral type WC8 + 07.5 

Marco fe S chmutz 1999) whose orbit is well determined, 
with a period of 78.53 ± 0.01 days, e = 0.326 ± 0.01, 
y W R = 68° ± 4° iSchmutz et al J 1 19971) and i = 63° ± 8° 
l|De Marco fe Schmutzlll999lh 

In a system consisting of a WR star with an early-type 
companion, a substantial contribution to the X-ray emission 
should come fro m the collision of the s tars' dense, highly su- 
personic winds <ICherepashchuklll976h . However, the situa- 
tion is complicated by the fact that the individual stars may 
also have substantial intrinsic X-ray emission from shocks 
that develo p in the winds as a result of lin e-driven insta- 
bilities fe.g. lOwocki. Castor fe Rvbickill 988h An important 
signature of colliding wind emission is phase-locked variabil- 
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ity of the flux and/or hardness of the X-ray emission. The 
simplest explanation for the variability of these observable 
quantities is the variation in the amount of absorption the X- 
rays suffer as the system moves through its orbit. Increasing 
absorption reduces the flux, but it also increases the hard- 
ness because softer X-rays are more strongly absorbed. In 
an eccentric binary, there may be other factors which give 
rise to phase- locked variability. Near periastron, the intrinsic 
X-ray luminosity may increase because the winds are denser 
when they collide. The emission will also be softer if the 
winds are not travelling at their terminal velocities, as the 
shock speeds will be slower near periastron compared with 
apastron. 

Einstein observations of ■y 2 Vel found that it was not 
unusually brig ht in X-rays, with ix/iboi(0) — 0.44 x 10~ 7 
iPollocklll98 7) [cf. the mean Lx/iboi for single stars in the 
Einstei n catalogue o f O stars is 2.5 x 10" " 7 JChlebowskll989l: 
see also[Moffat et all2002ft ]. There was however a suggestion 
that the X-rays from 7 2 Vel were soft er when the O star is 
in front of the WR star <Pollocldll987h . 

The first convincing evidence of the wind- wind collision 
in 7 2 Vel came not from X-ray data, but from International 
Ultraviolet Explorer (IUE) and Copernicus ultraviolet spec- 
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tra lISt-Louis. Willis fc Stevenslll993l) . They found that in 
broad terms the variability of the UV line profiles could be 
explained in terms of selective line eclipses of the star light 
by the WR star wind, with the effect being restricted to res- 
onance and low-excitation transitions of species common in 
the WC8 wind (e.g. Cll, Cm, Civ, Sim, Siiv and Siv). 
However, the details of some of the variability could not be 
explained by a simple spherically symmetric WR star wind. 
Instead there was evidence from resonance lines expected to 
occur in the winds of both stars (C IV, Si IV and S iv) that 
the collision between the stars' winds forms a cavity in the 
WR star wind. At phases when the O star is in front, one is 
looking into the cavity and sees the O star wind in absorp- 
tion. As the O star wind is ~1000 km s _1 faster than the 
WR star wind, a high- velocity blue absorption wing appears 
in the P Cygni profile. At phases when the WR star is in 
front, the O star wind travelling towards the observer is pre- 
vented from reaching its terminal velocity by the wind-wind 
collision, and so the high-velocity wing disappears, masked 
by the absorption trough of the WR star wind. 

X-ray emission from the wind-wind collision in 7 2 Vel 
was first studied in detail with ROSAT. F r om an analysis of 
13 observations, IWillis. Schild fc Stevensl (Il995l) found sig- 
nificant phase-dependent X-ray variability that is repeatable 
with binary phase. They found two main components to the 
emission. When the O star was not in front of the WR star, 
the emission in the 0.1-2.5 keV range was relatively constant 
and soft (kT ~ 0.19 keV) with L x ~ 2.5 x 10 31 erg s" 1 . 
When the O star was in front they found the X-ray emission 
was enhanced by the addition of a harder component, with 
kT> 1-2 keV and L x > 10 32 erg s" 1 . They attributed this 
harder component to emission from the wind-wind collision, 
observed through the cavity in the WR star wind formed 
by the O star wind. This was supported by hydrodynamical 
modelling of the wind-wind collision, which found that at 
ROSAT energies the emission from the wind- wind collision 
is only observable at phases when the O star is in front. 
Using the winds' terminal velocities in the hydrodynami- 
cal simulations overpredicted the luminosity in the ROSAT 
band (0.1-2.5 keV) by about an order of magnitude. Better 
agreement with the observations was obtained by assuming 
the O star wind collided at less than its terminal velocity, 
which is expected as the wind-wind collision region is likely 
to be close to the O star (as its wind is much weaker than 
that of the WR star). 

The hydrodynamical modelling in I Willis et alJ (Il995l) 
predicted copious X-ray emission at energies greater than 
2.5 keV , inaccessible b y RO SAT, but within the range of 
A SO A . IStevens et all Il99fJ) obtained two A SCA spectra 
taken near periastron (at phases $ = 0.978 and $ = 0.078, 
where $ = corresponds to periastron). These were com- 
pared with a set of synthetic X-ray spectra, generated from 
hydrodynamical simulations with a range of different wind 
parameters (namely the mass-loss rates and terminal ve- 
locities of the two winds). B y fitting their synth etic spec- 
tra to their observed spectra, Istevens et alJ (^996j) we re not 
only able to confirm the conclusion of I Wffh^etall (Il995l) 
that 7 2 Vel is a colliding wind system, but were also able 
to put constraints on some of the stars' wind parameters. 
Most notably, they found the mass-loss rate of the W-R 
star to be ~ 3 x 10 -5 Mq yr _1 . This is a factor of three 
lower than the mass-loss rate derived from radio observa- 



tions (8.8 x 10~ 5 M B vr -1 lBarlow. Roche fc Aitkenlll983) . 
IStevens et alJ lll996l) suggest this may be because radio ob- 
servations tend to overestimate mass-loss rates if the wind 
is i nhomogeneou s . 

iRauw et alJ i2000l) analysed an additional ASCA spec- 
trum of 7 2 Vel, this time ne ar apastron (< £ = 0. 570), and 
also reanalysed the data in IStevens et alJ lll996l). Having 
taken into account background sources. iRauw^tall (120001) 
effectively fitted the hard variable emission with a single 
temperature model. The bulk of the variability can be at- 
tributed to a changing column density towards the source. 
However, the temperature and luminosity also vary. The 
$ = 0.978 spectrum has a lower temperature (kT ~ 1.9 keV, 
using non-solar abundances) than the $ = 0.570 spectrum 
(kT ~ 2.7 keV). This could be because near apastron the 
winds are moving more quickly when they collide, result- 
ing in a larger post-shock temperature. The softest (and 
intrinsically most luminous) spectrum is seen at $ = 0.078 
(kT ~ 1.3 keV). If the soft X-rays are from the wind- wind 
collision, then at this phase they are being seen through the 
cavity in the WR star wind caused by the wind collision 
(whose position has been deflected by the Coriolis force). 
Alternatively, it could indicate that the intrinsic emission 
from the O star wind is making a significant contribution. 
In either case, at the other two phases the soft emission is 
absorbed by the wind of the WR star. 

The unprecedented spectral resolution offered by the 
Chandra Low- and High-Energy Transmission Grating Spec- 
trometers (LETGS and HETGS) and the XMM-Newton Re- 
flection Grating Spectrometer (E/AE ~ 100-1000) give us 
the opportunity to measure X-ray line shifts and widths to 
accuracies down to a few hundred km s _1 . This enables us 
to probe in detail the dynamics of the X-ray-emitting plasma 
in colliding wind binaries, giving new insights into the struc- 
ture of the wind- wind interaction regions in such systems. 

We present a new analysis of an archived Chandra 
HETGS ob servation of 7 2 Vel. T hese data have already been 
analysed bv lSkinner et alJ (1200 ll) . They found that the X-ray 
emission lines were broadened [full width at half-maximum 
(FWHM) ~ 1000 km s _1 ] but unshifted from their rest 
wavelengths. Furthermore, their results imply that the Ne IX 
line emission originates tens of stellar radii from the O star, 
well away from the central wind-wind collision region near 
the line of centres. Our analysis i ncludes a large r numb er of 
emission lines than published bv lSkinner et alJ 1120011) . We 
also discuss in more detail the interpretation of our results in 
terms of a colliding wind picture. The details of the observa- 
tion and the data reduction are described in Section [5] The 
broadband spectral properties are described in Section |3] 
while the results of fitting to the individual emission lines 
are described in Section 2] In Section we describe our at- 
tempts to model the emission line profiles using the model of 
iHenlev. Stevens fc Pittardl i2003T) . Our results are discussed 
in Section [(J and summarized in Section [7| Throughout this 
paper, the quoted errors are la confidence ranges. 



2 OBSERVATION AND DATA REDUCTION 

7 2 Vel was observed by the Chandra HETGS between 2000 
March 15 09:20:04 UT and 2000 March 16 04:08:12 UT, giv- 
ing 64 848 s of useful exposure time. The observation covered 
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orbital phases $ = 0.080-0.090, where <E> = corresponds 
to per iastron passage, using the ephemeris of lSchmutz et alJ 
(1997): 

JD (periastron) = 2 450 120.5 + 78.53-E (1) 

where E is the number of orbits since the periastron pas- 
sage on JD 2 450 120.5. Mid-observation was 6.7 days after 
periastron, and about 4 days after the passage of the O star 
in front of the WR star ($ = 0.03). The separation of the 
stars during the observation was 0.9 2-0.94 A.U. (assuming 
M W r = 9.5M(7) and Mp = 30M Q ; iDe Marco fc Schmut J 
ll999l:lDe Marco et aljEoOpl . 

The dispersed photons from the Chandra HETGS grat- 
ing assembly form a shallow X in the focal plane (see 
Fig. 0. One 'arm' of the X is from the High-Energy Grat- 
ing (HEG), the other is from the Medium- Energy Grating 
(MEG). The HEG has higher spectral resolution than the 
MEG: AA = 12 mA FWHM over the range 1.2-15 A for the 
HEG, versus AA = 23 mA FWHM over the range 2.5-31 A 
for the MEG {Chandra Proposers' Observatory Guide , Sec- 
tion 8.1). On the other hand, the MEG gives a higher signal- 
to-noise ratio in the first-order spectra of j 2 Vel. 

The data were reduced from the Level 1 events file us- 
ing CIAO v3.1, following threads available from the Chandra 
website 2 . This was done in preference to using the spectrum 
file obtained from the Chandra archive because it meant 
a more up-to-date version of CALDB (v2.27) could be ap- 
plied to the data. Similarly, response files (ARFs and RMFs) 
appropriate for this observation were generated following 
Chandra threads, in preference to using 'off-the-shelf re- 
sponse files. 

There are a total of 7850 counts in the non-background- 
subtracted first-order HEG spectrum and 16003 counts in 
the non-background-subtracted first-order MEG spectrum. 
For the analysis described here, the data were not back- 
ground subtracted, nor was the background spectrum mod- 
elled separately. For each of the two HETGS gratings, the 
background counts are extracted from two regions either side 
of and 4.5 times as wide as the source extraction regions 
(which are 4.8 arcsec wide in the cross-dispersion direction). 
There is therefore a possibility that photons from bright 
spectral lines could spread out into the background extrac- 
tion regions (because of the point-spread function of the 
telescope). Indeed, there is some evidence that the brighter 
emission lines in the source spectrum are also in the back- 
ground spectrum. This means that background subtrac- 
tion could adversely affect the measured results. Further- 
more, the Cash statistic, which was used for the analysis 
of emission lines described in Section 2] cannot be used with 
background-subtracted data. However, the count rates in the 
background extraction regions are less than 1 per cent of the 
rates in the corresponding source extraction regions (taking 
into account the fact that the background counts come from 
a detector area 9 times larger than the source counts). 



1 http:/ /cxc. harvard.edu/proposcr/POG/indcx. html 

2 http:/ /cxc. harvard.edu/ciao/threads/gspec. html 



3 BROAD-BAND SPECTRAL PROPERTIES 

The HEG and MEG spectra of 7 2 Vel are shown in Fig. H 
For illustrative purposes, the +1 and —1 orders have been 
co-added, and the spectra have been binned up to 0.02 A. 
The spectrum of j 2 Vel is dominated by strong emission lines 
from S (A « 4-5 A), Si (A « 5-7 A) and Mg (A » 7-9 A), 
with weaker lines from Ne (A ~ 9-14 A) and Fe (e.g. A ~ 2 A 
and A w 10 A). In this section we describe the analysis of 
the global properties of the spectrum, while in the following 
section we describe our analysis of the individual emission 
lines. 

The broad-band spectral analysis was carried out using 
XSPEC 3 vll.3.1. For the purposes of this analysis, the +1 
and —1 orders of each spectrum were co- added, and the 
data were binned so there were at least 20 counts per bin. 
This meant that the \ 2 statistic could be used, and the 
goodness of various spectral models could be assessed. The 
spectral models were fitted to the HEG and MEG spectra 
simultaneously, using all the available data. 

The ARFs used in the broad-band spectral fitting were 
generated by co-adding the ARFs for the relevant +1 and 

— 1 orders. A similar procedure cannot be followed for the 
RMFs, and so one must use the RMF for either the +1 or 

— 1 order for the whole co-added spectrum. Experimenta- 
tion showed that the combination of HEG and MEG RMFs 
used did not have a significant effect on the results. How- 
ever, it was noted that the value of \ 2 obtained when us- 
ing the HEG -1 RMF and the MEG +1 RMF was always 
lower than when using any of the other three possible com- 
binations. The differences in \ 2 were always small, but it 
may be connected to the fact that the HEG —1 order and 
the MEG +1 order contain the larger number of counts in 
their respective spectra [3951 counts (HEG —1) versus 3899 
counts (HEG +1) and 8543 counts (MEG +1) versus 7460 
counts (MEG — 1)]. The results described hereafter were ob- 
tained using this pair of RMFs. 



3.1 2T models 

A visual inspection of the Chandra HETGS spectrum of 
7 2 Vel indicates there is a wide range of temperatures 
present in the X-ray-emitting plasma, as evidenced by emis- 
sion lines from a wide range of ions [from Neix, whose 
emission peaks at T max m 4 MK (fcT max = 0.34 keV), to 
S xvi, whose emission peaks at T max « 25 MK (fcT max = 

2.2 keV)]. Hydrodynamical simulations of the wind- wind 
collision also predict a wide range of temperatures (see 
Section [SJ. One would therefore not expect a 2T thermal 
plasma model to provide a good fit to the spectrum; at 
best it can only characterize the emission, rather than ac- 
curately describe in detail the temperature structure of the 
X-ray-emitting gas. Furthermore, because of the wealth of 
detail in a grating spectrum (as opposed to a CCD spec- 
trum), a plasma emission model may have difficulty ac- 
curately fitting all the emission lines. The problems that 
arise could be due to several factors: Doppler shifting of 
lines, line broadening (on top of thermal Doppler broad- 
ening and instrumental broadening), non-solar abundances, 

3 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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Figure 1. The Chandra ACIS-S detector showing the dispersed photons of the HETGS spectrum of ■y 2 Vel. One arm of the X is from 
the HEG, the other is from the MEG. For the purposes of this figure, the data have been binned up by a factor of 16 to exaggerate the 
grating arms, and the image has been rotated so the ACIS-S detector is horizontal. 



and inaccurate atomic physics parameters (such as tran- 
sition rates and lab wavelengths) in the plasma emission 
model. Nevertheless, with these caveats in mind, fitting a 
2T model to the spectrum of j 2 Vel can give some insight 
into the general properties of the X-ray-emitting plasma. 
For this purpose, we used the xspec apec model (Astro- 
physical Plasma E mission Code 4 ; Smith & Bric khousel200fj; 
ISmith et al.ll200ll) . To model the absorption by the stellar 
winds we us ed the xspec wabs model (wh ich uses cross- 
sections from lMorrison fe McCammonl ll983'). 

We expect a significant portion of the X-ray emission 
to originate from the wind-wind interaction region, which is 
an approximately conical region around the O star whose 
apex points along the line of centres towards the WR star. 
The hot shocked gas streams along this cone away from 
the apex. At the time of the Chandra observation, the O 
star was approximately in front of the WR star, and hence 
the opening of the shock cone was towards the observer. 
There was therefore hot gas streaming along the shock cone 
towards the observer, and so we would expect to observe 
blueshifted emission lines. Because the opening of the shock 
cone was not directly towards the observer, there would in 
fact have been a range of line-of-sight velocities (distributed 
azimuthally around the shock cone), and so we would also 
expect to observe broadened emission lines. This simple pic- 
ture of the X-ray emission from the wind-wind collision is 
developed further in Section [4.51 

The amount of Doppler shifting and Doppler broaden- 
ing of the X-ray emission lines can be determined to a cer- 
tain degree from the broad-band spectrum (though this is 
done in more detail by investigating the individual lines; see 
Section 0. Line Doppler shifts can be modelled by thaw- 
ing the apec model's redshift parameter z. Note that this 
parameter measures the line-of-sight velocity of the X-ray 
emitting plasma, not the systemic line-of-sight velocity (as 
the shocked gas is moving with respect to the stars). Line 
broadening can be modelled by using the XSPEC gsmooth 
model, which convolves the whole spectrum with a Gaus- 
sian. The standard deviation a of this Gaussian can vary as 
a power-law with photon energy E: 

where Of, is the standard deviation at 6 keV. However, in 
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practice it was found that thawing the index a did not sig- 
nificantly improve the fit, so it was frozen at 1 (which gives a 
constant line width expressed as a velocity). For the fitting, 
one gsmooth component was used for the whole spectrum, 
and also the redshifts of the two apec models were con- 
strained to be equal (so there is one redshift characterizing 
the whole spectrum). 

We used non-solar abundances in our fitting, using 
the XSPEC vapec model. This model calculates abundances 
relative to hydrogen, which is essentially absent in WR 
winds. We overcame this by setting a large helium abun- 
dance [AT(He)/-V(H) = 10 6 , where N(X) denotes the num- 
ber abundance of element X], and then measuring abun- 
dances relative to helium. The abund ances of C, N and 
Q were fixed at _V(C)/_V(H e) = 0.14 .Morris et al]ll999_ 
ISchmutz fe De Marcel ll99Sl) . N(N)/N(He) = 1.0 x 10" 4 
■Llovd fe Stickland|ll999irand N(Q)/N( Re) = 0.028 [from 
N(C)/N(Q) = 5: iDe Marco et alJl2000| . Furthermore, as 
a starting point we used _V(Ne)/iV(He) = 3 .5 x 10" 3 and 
JV(S)/-V(He) = 6 x 10" 5 .Dessart et alJl200fl . and we fixed 
the relative abundan ces of all elements heavie r than neon 
at their solar values jAnders fe Grevessdll989l) . From this 
starting point, we allowed the abundances of Ne, Mg, Si, S 
and Fe to vary. 

The results of fitting the Chandra HETGS spec- 
trum of 7 2 Vel with an absorbed two-temperature 
thermal plasma model with non-solar abundances 
[wabs*gsmooth* (vapec+vapec)] are shown in Table 
Model A has a® and z frozen at 0, whereas in model B they 
are both free parameters. The data and the best-fitting 
models (concentrating on the S and Si line region) are 
shown in Figs. |3] (model A) and 01 (model B). Thawing 
OS and z significantly improves \ 2 1 although the fit is 
still formally unacceptable (the reduced \ 2 is 1-78 for 917 
degrees of freedom) . Figs. |3] and 0] indicate that the poor 
fits are mainly due to the emission lines, rather than the 
continuum. Nevertheless, the derived parameters can give us 
an overview of the general properties of the X-ray-emitting 
plasma. The 2 temperatures derived from the fit («9 and 
25 MK) are sensible given the range of peak emission tem- 
peratures of the lines observed in the spectrum. The value of 
ere corresponds to a FWHM of 1250^ km s _1 . The small 
measured line shift (z = —18 km s _1 ) implies the lines 
are not systematically shifted from their rest wavelengths. 
The values of eg and z can be compared with the results 
of fitting to individual emission lines (see Section 14.211 . 
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Figure 2. The first-order HEG and MEG spectra of -y 2 
up to 0.02 A. 



Velorum. The +1 and —1 orders have been co-added, and the spectra binned 



Thawing ctb and z does not affect the temperatures and 
emission measures of the two components. Also, while it 
seems to significantly affect the column density Nh, this 
is probably not the case because the errors in Table arc 
likely to be underestimated (given the poorness of the fits). 

While the abundances for a given model in TableQseem 
well constrained, a comparison between models indicates 
that this is in fact not the case. Furthermore, some mea- 
surements are hampered by the fact that for some elements 
only a few lines are available to constrain the abundance. 
For example, the neon abundance is really only being deter- 



mined from the Ne x Lya and Ly/3 lines, and this will lead to 
a temperature-abundance degeneracy. Indeed, uncertainties 
in the temperature distribution will affect all the abundance 
determinations. Unfortunately, this therefore means that the 
abundances in Table are rather unreliable. 

Note that the X-ray luminosities in Table Q have 
been deriv ed using the ■ffipjoarc os-determined distance of 
258+3^ pc iSchaerer et al.l fl99^) . This result has recently 
been thrown into doubt by the serendipitous discovery of 
an association of low-mass pre-main sequence (PMS) stars 



in the direction of 7 Vel IIPozzo et al 



pre- main sequence (FMb) stars 
(IPozzo et al J hoOCh . IPozzo et alJ 
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Table 1. Results of fitting the Chandra HETGS spectrum of ■y 2 Vel with an absorbed two-temperature thermal plasma model 
[models A and B: wabs*gsmooth* (vapec+vapec)] and an absorbed differential emission measure thermal plasma model [model C: 
wabs*gsmooth*c6pvmkl (see Section 13.21 1. Model A has [see equation and the redshift frozen at 0; in models B and C they 
are free parameters. The emission measures are expressed in terms of the helium number density rajjc (EM = J n a nu P dV), because of the 
lack o f hydrogen in the WR star wind. The quoted abundances are number ratios, relative to the solar number ratios lAnders fc Qrevessd 
1989). The X-ray luminosities are for the 0.4-10 keV energy range. 



Model ABC 



N n (10 22 cm" 2 ) 


1.92 ± 0.04 


2.08 ± 0.03 


2 24+ - 07 
-0.04 


fcTi (keV) 


n 77-3+0.012 
U - ' '''-0.014 


763+ 010 
u. /oo_ Q 0Qg 


EMi (10 54 cm' 3 ) 


2.29 ±0.08 


^«_0.15 




kT 2 (keV) 


2 19+ 008 


9 i c+0.04 
Z ' iO -0.05 




EM 2 (10 54 cm" 3 ) 


1.07 ±0.05 


96+ 07 

u - yu _0.03 




<T6 (eV) 


(frozen) 


10.6t°-2 


n O + 0.3 
°-°-0A 


Redshift z (km s ) 


(frozen) 


-18.3+H 


-160±i° 


(Ne/He) / (Ne/He) 


0.98 ± 0.12 


i S1 +0.19 


i ro+0.17 
J ~ oa -0.16 


(Mg/He) / (Mg/He) 


0.53 ± 0.03 


1 - ua -0.04 


i 1O+0.04 
1 ' 1J -0.06 


(Si/He) / (Si/He) Q 


n oc+0.02 
U ' 8t, -0.03 


1.34 ±0.03 


1 24+ 03 


(S/He) / (S/He) Q 


1.17 ±0.07 


1 66+ 07 


1.50 ±0.07 


(Fe/He) / (Fe/He) 


0.55 ± 0.06 


og+0.06 


n fis+ 06 

U - b8 -0.05 


Xl (d.o.f.) 


2.81 (919) 


1.78 (917) 


1.95 (914) 


L abs ( 1Q 32 erg g -l) 


1.0 


1.1 


1.1 


L£* (10 32 erg s" 1 ) 


6.7 


8.4 





( 2000) argue that these PMS stars are at approximately 
the same distance and age as y 2 Vel, placing 7 2 Vel within 
the Vela OB2 association at a distance of 360-490 pc, in 
good agreement with older distance estimates (e.g. 460 pc; 
IConti fc Smithlll972l) . If this larger distance is correct, the 
X-ray luminosities in TableQshould be increased by a factor 
of ~3. 

The absorbing column Ah derived from the fit corre- 
sponds to a vi sual extinct i on A v « 9, using the empir- 
ical relation of iGorensteinl dl975l) . In contrast, the visual 
extinction derived fr om optical studies is A v = 0.00-0.12 
l)van der HuchdfeoOll ). This implies that there is extra ab- 
sorption local to the source, in addition to the interstellar ab- 
sorption. Given the orientation of the orbit and the phase of 
the observation, the wind-wind collision zone was observed 
through the wind of the O star. The value of Ah is consistent 
with the column density through the wind of the O star to a 
point in the wind-wind collision zone that lies on the line of 
centres, which may reasonably be assumed to be a charac- 
teristic column density for the whole emission region. This 
can be shown by considering the optical depth t„ through a 
spherically symmetric wind expanding at constant velocity 
Woo to a point at (r,9,(f>) (in spherical coordinates centred on 
the star, with the z-a xis along the line of sight), which is 
given by dlgnacell200ll) 



hydrogen and mn is the mass of a hydrogen atom. Thus, by 
combining equations @ and we obtain 



(3) 



where _R* is the stellar radius and 7y = k v M /AnVpo R* is a 
characteristic wind optical depth JOwocki fc Cohenll200ll) . 
where M is the stellar mass-loss rate and k„ is the opacity. 
The optical depth is related to the column density Ah by 



t v = J K^pdr ps / K^munndr — K^muNn (4) 
where p is the mass density, tih is the number density of 



AT H = 



M 



AirmytVoor sin 9 



(5) 



The angle 9 [that is, the polar angle (in this coordinate sys- 
tem) of a point on the line of centres] is related to the orbital 
phase angle 9 (ty = 0° corresponding to the O star being in 
front) and the orbital inclination i by 

cos 9 = — cos 9 sin i (6) 

[See also equation (4) in iHenlev et alJ i2003t) : note that 
the 9 used in that equation is the supplement of the 9 
in equatio ns (0 and JSJ.] Using the orbital elements of 
ISchmutz et all (Il997t) (e = 0.326, wwr = 68°) and the phase 
of the Chandra observation ($ = 0.085 at mid-observation; 
see Section |21 , we obtain ^ = 36° . Combining this with 
i = 63° JPeMarco fc Schmutdll999l) gives 6 = 136°. Us- 
ing the separation (~0.93 A.U. ) and the wind mome ntum 
ratio (MwR"WR/Afo"o ~ 33; iDe Marco et al.ll2000l) . the 
distance along the line of centres from the O star to the 
wind-wind collision zone is r ~ 2 x 10 12 cm. The mass- 
loss rate and wind velocity are more difficult to determine, 
and the situation is complicated by the fact that the wind 
will not be at its terminal velocity at the wind collision. 
One consequence of this is that the column density will be 
larger than that given by equation ijjjj, because the inner re- 
gions of the wind are denser than is assumed by a constant- 
expansion model. Assuming M = 0.1-1 x 10~ 6 M yr _1 and 
Woo = 2000 km s _ (which are sensible values for a late O 
star), we find Ah ~ 0.3-3 x 10 22 cm -2 , which is consistent 
with the values derived from the spectral fitting. However, it 
should be noted that the wabs absorption model in xspec as- 
sumes neutral absorbing material. This will overestimate the 
opacity of the partially ionized stellar wind, and so underes- 
timate Nil- This discr epancy will be wors t at lower energies 
( < 1 keV; see fig. 4 in lCohen et"allll996l) . 
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Figure 3. The binned HEG (top) and MEG (bottom) spectra from 7 2 Vel plotted with the best-fitting absorbed 2T apec model, with 
the redshift z and erg frozen at (model A in Table 0. The residuals are the direct differences between the data and the model (in 
normalized counts s — 1 A -1 ). 



If we invert the above argument, and specify our mea- 
sured value of Nu in equation we obtain Mo ~ 8 x 
10~ 7 Mq yr~ . However, this value may not be very reli- 
able, given the approximations used to obtain it. 



3.2 DEM modelling 

As has already been stated, a 2T apec model can only char- 
acterize the spectrum, rather than accurately describing the 
temperature structure of the X-ray-emitting plasma. One 
would expect to obtain a more detailed picture of the tem- 



perature structure by using a differential emission measure 
(DEM) model. 



ISkinner etal] fcOOll) fitted the first-order MEG spec- 
trum of 7 2 Vel with an absorbed differential emission mea- 
sure (DEM) model, and found a strong emission measure 
peak near ~9-10 MK. We attempted to reproduce this, us- 
ing a wabs*gsmooth*c6pvmkl model in xspec. This model 
is based on the mekal thermal plasma model (unfortunately 
there is no DEM model in xspec based on apec). The DEM 
is expressed as an exponential of a polynomial 
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Figure 4. As Fig. [3] but with the redshift z and C6 free to vary (model B in Table [I]. 



DEM(T) oc exp 



^> fc P fc (T) 



(7) 



where Pk is a Chebyshev polynomial of order k, and the co- 
efficients afe are free to vary during the fit. As in the previous 
section, the abundances of Ne, Mg, Si, S and Fe were free 
parameters, as were oe and z. Again, we fitted the HEG and 
MEG spectra simultaneously. 

Our best- fitting DEM is shown in Fig.|5] and the result- 
ing spectral parameters are in Tabled (model C). The DEM 
mekal model actually gives a worse fit than the 2T apec 
model. This is most likely due to the fact that mekal is an 
older code, and so may have inaccurate wavelengths, tran- 
sition rates, etc. compared with the more up-to-date apec 



model. This means the derived spectral parameters should 
be treated with caution (in particular the redshift z, which 
is significantly different from that obtained in Section f3.H . 

Our DEM is broadly peaked at ~7-8 MK, compared 
with t he strong peak at ~10 MK found by ISkinner et all 
teOOlf) . To estimate the uncertainty on the shape of the 
DEM, we have calculated the errors on the coefficients 
using XSPEC, and then used these in a Monte-Carlo simula- 
tion to generate a set of 1000 DEMs. While the small peak 
at ~35 MK does not seem to be very significant, the drops 
in the DEM below 7 MK and above 40 MK do seem to be 
real. The DEM therefore indicates a wide range of plasma 
temperatures from a few MK to ~40 MK, as expected from 
the range of emission lines seen in the spectrum. Without 
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Figure 5. Our best fit differential emission measure (DEM) 
model for 7 2 Vel. The arrows (at 8.85 MK and 24.9 MK) indicate 
the best fit temperatures from our 2T apec model fitting (TabletTl 
model B). The rise in the DEM below ~2 MK is probably an arte- 
fact of the assumed functional form of the DEM [equation 1711 . 

knowing more about ISkinner et alf s j200lTl DEM, we can- 
not say whether or not our DEM is consistent with theirs. 

An alternative method of calculating the DEM (with 
uncertainties) would be to use the method described by 
IWoidowski fc Schubd j2004t) , which involves fitting the spec- 
trum with a number of components, each of which is the 
entire X-ray emission line spectrum for a single ion. The 
DEM is thus pieced together from the emission measures of 
these individual components (each of which corresponds to 
a different temperature range). However, such an additional 
analysis of the DEM is beyond the scope of this paper. 

Rather surprisingly, our DEM also shows a pronounced 
rise below ~2 MK. However, as there is very little emis- 
sion above ~14 A, the cool end of the DEM (below ~4 MK, 
the peak emission temperature of Ne ix) will be very poorly 
constrained. This rise in the DEM is thus almost certainly 
just an artefact of the assumed functional form of the DEM 
[equation 10]. An unfortunate consequence of the rise in the 
DEM at low temperatures is that the intrinsic (unabsorbed) 
flux is over-estimated by several orders of magnitude at low 
energies. For this reason we have not quoted a value for Lx 
for this model in Tabled This excessive intrinsic low-energy 
flux probably explains why the column density derived from 
the DEM model is somewhat higher than that derived from 
the 2T model (see Table [TJ. 



4 EMISSION LINE PROPERTIES 

Unlike the broad-band spectral analysis, the analysis of the 
individual lines in the Chandra spectrum of j 2 Vel was car- 
ried out using unbinned, non-co-added spectra, so no spec- 
tral information was lost. Because so me bins con tained low 
number of counts, the Cash statistic llCashlll979l) was used 
instead of x 2 - This analysis was carried out using SHERPA 
(as distributed with CIAO v3.1). 

The analysis of the X-ray spectrum concentrated on 
the emission lines from H- and He-like ions of Ne, Mg, 
Si and S. For all four of these elements the H-like Lya 
lines and He-like fir (forbidden-intercombination-resonance) 
triplets were detected. Ly/3 lines were also analysed for 



Nex, Mgxn and Sixiv, but not Sxvi, as its Ly/3 line is 
blended with the Arxvn fir triplet near 3.95 A. The Nex 
Lyman series extends to include Ly7 and Ly<5. He/3 lines 
(ls3p 'Pi -> Is 2 %) were also detected from Sixni and 
Mg XI. 

Each emission line or multiplet was analysed individ- 
ually over a narrow range of wavelengths (typically 0.2- 
0.5 A). The fitted model consisted of a continuum compo- 
nent plus Gaussian component (s) to model the line emission. 
Although colliding wind binaries are expected to e xhibit a 
wide range of line profile shapes (iHenlev et alj l2003). Gaus- 
sians give good fits to the emission lines in j 2 Vel (see Sec- 
tion 14. H , and thus provide a good way of quantifying the 
line shifts and widths. It was found that the choice of con- 
tinuum model (constant or power-law) did not significantly 
affect the best-fit parameters of line emission component. As 
a result the simpler constant model was chosen in preference. 
The nature of the line emission component used depended 
on the nature of the line being investigated, as discussed 
below. 

Lyman lines are in fact closely spaced doublets due to 
the spin-orbit splitting of the upper level. The splitting is 
~5 mA for Lya and ~1 mA for Ly/3. This separation is too 
small to be resolved by the Chandra gratings, and so a good 
fit can be obtained with a single Gaussian. However, it is 
not obvious what should be used as the rest wavelength of 
the line if one wishes to measure a Doppler shift. The Ly- 
man lines were therefore fitted with two Gaussians whose 
parameters were tied together. The Doppler shifts of the 
two components (expressed as velocities) were constrained 
to be equal, as were their widths. The relative intensities 
of the two components were fixed at the theoretical values 
from ATOMDB vl.1.0, the shorter-wavelength component be- 
ing approximately twice as bright as the longer-wavelength 
component. 

In contrast to the Lyman lines, the components of the 
He-like fir triplets are resolvable by the Chandra gratings. 
These lines were fitted with three Gaussians whose parame- 
ters were tied together (the intercombination line is in fact a 
closely spaced doublet, but using four Gaussians did not sig- 
nificantly alter the results). The Doppler shifts and widths 
were linked as for the Lyman lines. For this purpose, the lab- 
oratory wavelength of the brighter, longer-wavelength inter- 
combination line was used. The amplitudes of the Gaussians 
were all free to vary. The He/3 lines are singlets, and so these 
were fitted with single Gaussians. 

4.1 Fits to individual lines 

The fit results are shown in Table [5] All these results 
were obtained by fitting the line model to all four spectra 
(HEG —1 and +1, MEG —1 and +1) simultaneously, except 
for the Sxvi Lya line (the results were obtained just from 
the MEG data, as a sensible fit could not be obtained if 
the HEG data were included) and the Sixiv Lya line (the 
MEG —1 data were excluded from the fit; see below). For 
each individual spectrum the appropriate ARF and RMF 
was used. 

The Lya lines and He-like fir triplets from S, Si, Mg and 
Ne are shown in Fig. |5] along with the best-fitting model 
for each. For illustrative purposes only the data have been 
binned up to 0.01 A in these plots. For comparison, the un- 
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Table 2. The measured wavelengths (A bs)i widths (AA) and fluxes of the X-ray emission lines in the Chandra HETGS spectrum 
of 7 2 Vel. All results were obtained by fitting to the unbinned, non-co-added HEG and MEG spectra simultaneously (except where 
indicated). Parameters without quoted errors were tied to other fit parameters, rather than being free. Laboratory wavelengths (Ai a t,) 
and temperatures of maximum emission (T max ) are from ATOMDB vl.1.0. 
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a Results obtained just by fitting to MEG data 

b Results obtained by omitting MEG —1 data from the fit 

binned Sxv fir triplet and Sixiv Lya line are shown in 
Fig. |7| Note that there appears to be some excess emis- 
sion in the S xv fir HEG data approximately 0.01 A short- 
wards of the forbidden line. However, this is probably not 
real, because it does not appear at the same wavelength 
in the ±1 and —1 orders, it does not appear at all in the 
MEG data, and there is no corresponding feature short- 
wards of the resonance line. Note also that the MEG — 1 
Si xiv Lya line appears significantly different from the oth- 
ers - it is shifted towards the red, and also it appears 
to be skewed redwards. However, this skewing is proba- 
bly not real - if it were it should be detected in the other 
three spectra (if anything, the line in the MEG ±1 spec- 
trum appears to be slightly skewed bluewards). The wave- 
length obtained by fitting to just the MEG —1 spectrum 
is 6.1834 ± 0.0014 A, that obtained by fitting to the other 
three spectra is 6.1795 ± 0.0006 A, and that obtained by fit- 
ting to all four spectra is 6.1803±0.0006 A. Hence, while the 
MEG —1 spectrum gives a significantly different wavelength 
from the other three spectra, including it in the fit does not 
significantly affect the results (this is true of the FWHM and 



flux as well). Nevertheless, as the line is clearly different in 
the MEG —1 spectrum, we have chosen to quote the results 
obtained by omitting the MEG —1 spectrum from the fit in 
Table H 

The fit results are illustrated in Fig. |S| which shows 
the measured line shifts and widths as functions of labo- 
ratory wavelength. One can see from the results that the 
lines are typically unshifted from their lab wavelengths, with 
FWHMs of ~1000-1500 km s" 1 . The shifts and widths are 
uncorrelated with lab wavelength or the ionization potential 
of the emitting ion. The line shifts are well fit with a single 
mean shift (independent of wavelength) of — 64 ± 12 km s _1 
{xl = 1-13 for 14 de grees of freedom), while the line 
widths are adequately fit with a single mean FWHM of 
1240 ± 30 km s" 1 ( X l = 1.91 for 14 degrees of freedom). 

While the Cash statistic can be used on data with low 
numbers of counts, it does not contain any goodness-of-fit 
information. Therefore, to get some idea of the goodness of 
the fits, we again co-added the +1 and —1 orders of the HEG 
and MEG spectra, binned up the data so there were at least 
10 counts per bin, and then repeated the fitting procedure 
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Figure 6. The Lyct lines from Sxvi, Si XIV, Mgxn and Nex and the He-like fir triplets from Sxv, Sixin, Mgxi and Neix observed in 
the Chandra HETGS spectrum of 7 2 Vel. The solid line in each panel shows the data, while the dashed line shows the best fit model. 
For illustrative purposes only, the +1 and —1 orders of each grating (HEG and MEG) have been co-added and binned up to 0.01 A. 
Note that the Sxvi Lyct fit was just carried out on the MEG data, and so the HEG data are not plotted. 



using x 2 (which does enable us to assess the goodness of fit) . 
Just for this x 2 fitting, the HEG -1 and MEG +1 RMFs 
were used, as in the broad-band spectral fitting. To illustrate 
the effect of rebinning the data, the binned Sixiv Lya line 
is shown in Fig. El We found no significant difference in 
the wavelengths, widths and fluxes derived using the two 
different fit statistics - all agreed within la. In all cases, 
Gaussians give good fits to the emission lines. A possible 
exception to this is the Si xiii fir triplet - Gaussians do not 
give good fits to the HEG lines (which appear to be slightly 
skewed redwards) but they do to the MEG lines. However, 
this skewing is probably not real, as there is no evidence 
for it in the MEG data (which has higher signal-to-noise). 



There is therefore no convincing evidence that Gaussian line 
profiles do not give good fits to the observed emission lines. 

There is no improvement in \ 2 when two Gaussians are 
used for the Lyman lines as opposed to just one. However, 
while it is not required by the data, the use of a doublet 
model is justified on physical grounds, as discussed above. 



4.2 Comparison with broad-band results 

The results described in this section can be compared with 
the results derived from the broad-band spectrum in Sec- 
tion Fig. HOI compares the line shifts and widths derived 
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Figure 7. The unbinned Sxv fir triplet (left) and Si XIV Lya line (right), with the best fitting Gaussian line profiles. The anomalous 
appearance of the MEG —1 line in the right figure is discussed in the text. 
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Figure 8. Line shifts (left) and widths (right) as functions of laboratory wavelength. The vertical lines divide the plot up into lines from 
S, Si, Mg and Ne (from left to right). Lines from H-like and He-like ions are denoted with solid and dashed error bars, respectively. 



in this section with those inferred from the values of z and 
as in Table Q (model B). There is excellent agreement be- 
tween the width derived from the broadband fitting (FWHM 
= 1250+5° km s _1 ) and the mean width of the individual 
emission lines (FWHM = 1240 ± 40 km s" 1 ). Furthermore, 
the lack of any observed correlation between line width and 
wavelength justifies our use of a — 1 in equation J3J. 

We would expect there to be good agreement between 
the line shifts derived from the broadband fitting and the 
mean shift of the individual emission lines, as the rest 
wavelengths in the apec model used to measure the for- 
mer are the same as those used to calculate the latter. 
In fact, the line shift derived from the broadband fitting 
(z — — 18.3+J 5 km s -1 ) disagrees with the mean shift of 
the individual emission lines (—64 ± 12 km s _1 ) at the 3<r 
level. However, the agreement is probably better than this 
indicates, because the poor spectral fits in Section [3] mean 



the errors on the spectral parameters (including z) are likely 
to be underestimated. 

4.3 Estimating the temperature using line flux 
ratios 

There are two simple ways to estimate the temperature T 
of the X-ray-emitting plasma using the flux ratios of emis- 
sion lines. The first uses the flux ratio of the H-like Lya 
line to the He- like resonance line for a given element. As 
the temperature increases the ionization balance shifts from 
He- to H-like ions, and this change is manifested in the line 
flux ratio. The second uses the flux ratio G — (/ + i)/r of 
the resonance (r; ls2p x Pi — > Is 2 1 So), intercombination (i; 
ls2p 3 Pi, 2 -> Is 2 l S ) and forbidden (/; ls2s 3 Si -> Is 2 ^o) 
lines of a He- like ion. G is a decreasing function of tem - 
perature (e.g. iPorauet et alJliooll: iPaerels fc Kahnl 120031) . 
due to the different temperature dependencies of the mech- 
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Figure 10. Comparison between line shifts (left) and line widths (right) derived from fitting to the broad-band spectrum and those 
derived from fitting to individual emission lines. The data points in both figures are the results from fitting to individual emission lines. 
The dot-dashed line in the top figure shows the mean shift of the lines (—64 km s — 1 ). The solid line in the top figure shows the line shift 
corresponding to z = —18.3 km s" 1 , and the solid line in the bottom figure shows the line width corresponding to a@ = 10.6 eV (model 
B, Table |T}. 
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Figure 9. The Si XIV Lya line shown binned up so there are 
at least 10 counts per bin, with the best fit Gaussian line profile, 
obtained by fitting to the two binned spectra simultaneously using 



anisms for populating the upper levels of the lines (direct 
collisional excitation from the ground state for the reso- 
nance line, cascades from higher levels populated by dielec- 
tronic recombi nation for the forbidden and intercombination 
lines; see ISmith et alJl200llb The G ratio is also sensitive to 
electron density at high densities, due to population of the 
upper level of the r esonance line from the ls2s 1 So level 
llPorcmet et alJl200ll) . However, the d ensities at which this 
becomes important (n e > 10 12 cm -3 ; IPorauet et al.l l200ll) 
are larger than those expected in the wind-wind collision 
in 7 2 Vel (the hydrodynamical simulations in Section in- 
dicate ne^lO 11 cm -3 ). 

To estimate temperatures using H-like to He-like 



flux ratios, we use line emissivities calculated as func- 
tions of T with the Astrophys ical Plasma Emission Code 
and Database (a pec/aped 5 ; ISmith fc Brickhousel l200fJ : 
ISmith et alJ 1200 if) . These emissivities are distributed in 
ATOMDB vl.1.0. Comparing the measured H-like to He-like 
flux ratios to the calculated values, we infer temperatures of 
5.6 ± 0.5 MK for Ne, 8.6 ± 0.3 MK for Mg, 11.7 ± 0.2 MK 
for Si and 14 ± 1 MK for S. 

The measured G ratios are 0.9 ± 0.3, 0.67 ± 0.10 and 
0.91 ±0.05 for Neix, Mgxi and Sixm, respectively. To cal- 
culate tempe ratures, we compare t hese values to the values 
calculated bv lPorauet et all ^2001) in the low-density limit. 
Their calculations include the effects of blending of dielec- 
tronic satellite lines. The contribution of these blended lines 
to the measured fluxes depends on the spec tral resolution 
of the instrument used llPorauet et alJl200ll) . and so they 
tabulate G ratios as a function of T assuming the MEG 
spectral resolution [AA(FWHM) = 23 mA] and assuming 
the HEG spectral resolution (AA = 12 mA). Our fits use 
data from the MEG and the HEG simultaneously. However, 
this does not matter as there are no significant differences 
between the tempera tures inferred using th e MEG values or 
the HEG va lues from IPorauet et al] fcOOlT) . Using the MEG 
values from IPorauet et all (1200 ll)" we infer temperatures of 



3+3 



MK for Neix, 711 MK for M S XI and 4 -8-a7 MK for 
Si xiii, whereas using the HEG values we infer 2±% MK for 
Neix, 7± 3 2 MK for Mgxi and 5.4±°;® MK for Sixm. 

The temperatures measured using the fir triplets of 
Neix and Mgxi are in good agreement with those measured 
using the H-like to He-like flux ratios. However, the tem- 
perature measured from the Si xiii fir triplet is significantly 
cooler than that measured from the Sixiv:Sixill flux ratio. 
The calculation of the temperature from the G ratio may be 
affected by blending of the Si xiii forbidden line with Mg XII 
Ly7. This would mean the G ratio is over-estimated and T 
is under-estimated by this method, in agreement with the 
observed discrepancy. However, the discrepancy may more 
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simply be due to the He- like emission coming from different, 
cooler regions than the H-like emission. 

4.4 The location of the X-ray emitting plasma 

The flux ratio R = f/i is sensitive both to electron density 
n a and the UV radiation field, due to depopulation of the 
upper level of the for bidden line to the upper level of the in- 
tercombination line jGabriel fc Jordanlll969tlPorauet et alJ 
l200ll:lPaerels fc Kahnl200St . Comparing the measured R ra- 
tios with theoretical values gives information on the electron 
density and/or UV radiation field, and with some other as- 
sumptions the location of the X-ray emitting plasma can be 
inferred. 

The measured R ratios are 2.6 ± 1.5 f or Ne , 1.0 ± 0.2 
for Mg and 2.5 ± 0.2 for Si. iPorauet et alJ J200ll) have cal- 
culated R ratios as functions of n c for a range of electron 
temperatures, radiation te mperatures a nd radiation dilution 
factors [W(r) = 0.5(1 — yl — (R*/r) 2 ), where r is the dis- 
tance from the centre of a star of radius R,]. The radius of 
the WR star is 3.27? Q <De Marco et alJl2000h . and the dis- 
tance from the WR star to the wind- wind interaction will be 
at least half the separation (i.e. r > 1007?©). The dilution 
factor of the WR star's radiation field is thus W < 0.00026, 
and its effect on the observed R ratios will be negligible. 

Fig. II II shows the R ratio as a function of n D in the ab- 
sence of a radiation field (T ra d = 0), and in the presence of a 
blackbody radiation field w ith T ra d = 35 000 K (the effec tive 
temperature of the O star; iDe Marco fc Schmutdir999T) for 
two different values of the dilution factor, cal culated by av- 
eragin g the data for 30 000 K and 40 000 K in lPorauet et alJ 
(2001). The plotted curves are for electron temperatures 
of 3.0 MK (Neix), 6.3 MK (Mgxi) and 5.0 MK (Sixm). 
T hese temperatures we re chosen from the values tabulated 
bv lPorauet et all i200ll) based on the values derived in using 
G ratios in Section f4.3l though in fact the choice of electron 
temperature has no effect on the conclusions. Also shown in 
Fig. Illl are the measured R ratios. 

For Neix the observed value of the R ratio is con- 
sistent with the theoretical low-density, low-UV flux value 
(R « 3). Furthermore, we can see that even a relatively 
weak UV flux (W ~ 0.01) would give rise to an R ratio 
significantly lower than that which is observed. We there- 
fore conclu de that W <g 0.01. The r adius of the O star 
is 12.47? Q iDe Marco fc Schmutd ir999V and hence for the 
Ne ix-emitting plasma, r 2> 57?, (O) — Ax 10 12 cm. By the 
same argument, we can see from Fig. Illl that W < 0.01 for 
Si xiii, and so we have r > 57?, (O) for the Sixill-emitting 
plasma. 

For Mgxi the observed value of the R ratio is signifi- 
cantly lower than the theoretical low-density, low-UV flux 
value (R « 2.7). If the UV radiation field were negligi- 
ble, this would imply an electron density of ~10 13 cm -3 , 
which is far larger than is expected in the wind-wind colli- 
sion zone. Even with a dilution factor of W ~ 0.01, the im- 
plied density (~10 12 -10 13 cm 3 ) is too large. On the other 
hand, if the dilution factor were as large as W ~ 0.10, the 
value of the 7? ratio would be significantly lower than that 
which is observed. Hence for the Mg xi-emitting plasma, 
0.01 < W < 0.10, and so 1.77?»(0) < r < 5R*(0). This 
means the Mgxi emission emerges from closer to the O 
star than the Sixm emission. The implications of this are 



Ne IX 



4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




2.5 
2 

rr 1.5 

1 

0.5 



2.8 
2.6 
2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 



Trad = 

T rad = 35 kK, W = 0.01 
T raH = 35 kK,W = 0.10 



100 10000 



1 e+06 1 e+08 
n e (cm" 3 ) 

Mg XI 



1e+10 1e+12 1e+14 



Trad = 

T rad = 35 kK, W = 0.01 
T rad = 35 kK, W = 0.10 




100 10000 



1e+06 1e+08 
n e (cm" 3 ) 

Si XIII 



1e+10 1e+12 1e+14 













T rad = 






T rad - 35 kK, W - 0.01 






T rad = 35kK, W = 0.10 








- 



100 10000 1e+06 1e+08 1e+10 1e+12 1e+14 
n e (cm -3 ) 



Figure 11. Measured and calculated R ratios for Neix, Mgxi 
and Si XIII. The thin curves show the R ratio as a function of n e 
for different radiation temperatures (T ra[ i) and dilution factors 
(W). The thick horizontal lines show the measured R ratios and 
the la errors. 



discussed in Section |S| It should be noted, however, that 
the uncertainty in the O star's emergent photospheric flux 
at the far- and extreme-ultraviolet wavelengths relevant for 
the photoexcitation transiti ons (1270 A for Neix , 1033 A for 
Mgxi and 864 A for Si xill: lPorauet et alJl200lTl may be sig- 
nificant, and that a 35 000 K blackbody does not represent 
the flux of the O star very well, especially at the wavelengths 
relevant to the higher Z elements. A more accurate model of 
the photospheric flux will likely modify the curves in Fig. 1111 
in a way that is not easy to predict. 

ISkinner et alJ ||200 Jl use this dataset to derive distances 
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of the X-ray-emitting plasma from the O star by considering 
the following relationship between the R ratio, the electron 
density n c and the UV photoexcitation rate <j> from the / to 
i upper levels: 



1 + {4>/4> c ) + (n c /n c ) 

where n c is t he critical density and <j> c the critica l photoex- 
citation rate iBlumenthal. Drake fc Tuckerlll972fr . Ro is the 
low-density, low-UV flux l imit of the R r a tio, i. e. if n e <C n c 
and 4> "C <j>c, R — * Ro- ISkinner et alJ i200lT) argue that 
R > |7?o for their fir triplets (as lower values of the R ratio 
would be unambiguo usly detected), an d henc e <f>/4> c < 0.5. 
Using the results of iBlumenthal et alJ (I1972T) . they derive 
minimum line formation radii of r m i n ~ 3i?*(0), 9i?»(0) 
and 30ii, (O) for Sixm, Mgxi and Neix, respectively. Their 
results for Sixm and Neix are c onsistent with our res ults. 
Using t he more recent res u lts of IPorauet et alJ feOOll) in- 
stead of IBlumenthal et alJ (Il972l) . we have increased r m i n 
for Si xni from 3i?»(0) to 57?, (O). On the other hand, our 
analysis does n ot place a strong cons traint on r m i n for Neix, 
mainly because IPorauet et alJ <!200ll) do not quote values of 
the R ratio for dilution factors below 0.01. 

It is not surprising that Ne IX originates far from the O 
star. Near the line of centres (and hence near the O star) 
the winds collide head-on, producing shock-heated plasma 
that is too hot for Neix to exist (tens of MK). This plasma 
moves outwards along the shock-cone, cooling radiatively. It 
will not cool to the temperature at which Neix emission 
peaks («4 MK) until it has travelled some distance (the 
hydrodynamical simulations discussed in Section ^indicate 
the plasma cools to ~4 MK at a distance of ~10i?»(O) from 
the O star). This is a factor of 3 less than the mi nimum line 
formation radius derived bv lSkinner et all J2001) . However, 
an inspection of the models used in Section [3.11 shows that 
there is considerable iron emission near the Neix fir triplet; 
combined with the fact that there are very few counts in 
this part of the spectrum, this means the Neix results are 
probably untrustworthy. 

Our results for the locatio n of the Mgxi-em itting 
plasma are different from those of lSkinner et al.l (1200 ll) . We 
find this plas ma is lo c ated at 1.7R, (O) < r < 5R,(0), 
whereas Iskinner et"aH l|200l]) derive r > 9_R,(0). Part of 
this discrepancy may come from the fact we measure d iffcr- 
ent values of the R ratio for Mg xi: 1.8±j5;jj JSkinner et alJ 
1200 ll) versus 1.0 ± 0.2 (this work). This difference may be 
due to differences in the calibration used, or to the fact 
that lSkinner et all i200ll) use background-subtracted HEG 
data only, whereas we use data from both gratings without 
background subtraction. It should be noted that our mea- 
sured width for the Mgxi fir triplet is significantly larger 
than most of the other widths in Table |2] In order to as- 
sess whether or not this affects our flux measurements, we 
have repeated the fit with the triplet's FWHM fixed at 
1200 km s" 1 (the average FWHM of all the lines). These 
measured fluxes are ~ 10-20 per cent lower than our origi- 
nal measured values, but the differences are not significant. 
Furt hermore, the R ratio is unaffected. 

ISkinner et alJ i200ll) further state that uncertainties 
in their flux measurements and possible blending with 
Nex Lyman lines means that their value of the Mgxi R 
ratio could be consistent with the low-density limit Rq. 



The Nex Lyman lines from upper levels 6-10 have wave- 
lengths of 9.362, 9.291, 9.246, 9 .215 and 9.194 A, respec- 
tively llHuenemoerder et al.ll200lT) . whereas the wavelengths 
of the Mgxi forbidden and intercombination lines are 9.314 
and 9.231 A (atomdb vl.1.0). Given the wavelengths of 
these Ne x lines, and without any information on their emis- 
sivities (these lines are not in aped), it is difficult to say 
which of the / or i line will be more affected by blending, 
and hence one cannot say for sure that the measured value 
of the R ratio for Mgxi is probably underestimated. 

If we do assume that blending with Nex Lyman lines 
tends to lower the observed value of the R ratio for Mgxi, 
we can derive a location for the Mg xi-emitting plasma using 
the method in lSkinner et alJ <l200ll) . However, whereas they 
used R > |-Ro, given our results we can only say that R > 
|i?o- This implies 4>/4> c < 2.0, and gives r > 4.57?* (O), 
which is consistent w ith the range we derive above using the 
IPorauet etHl d200ll) data. 

4.5 Geometry of the wind-wind collision 

Some insight into the geometry of the wind-wind collision 
zone may be gained by assuming the X-ray line emission 
comes from a cone with opening half-angle whose symme- 
try axis lies along the lin e of centres with the apex point - 
ing towards the WR star <Liihrdll997t iPollock et alJl2004T) . 
This arrangement is illustrated in Fig. 1121 The angle 7 be- 
tween the line of centres and the line of sight is the sup- 
plement of the angle 6 in equation © and hence is given 
by cos 7 = cos \1/ sin i. We assume the emitting material is 
flowing along the cone at speed vo and that there is no az- 
i muthal velocity component [these are explicit assumptions 
of lLiihrJ dl997ll model]. The former is justified in the light 
of 2-D hydrodynamical simulations of the colliding winds 
(such as those in Section |SJ , which show that the velocity 
vector is indeed directed along the shock cone. There is a 
small divergence in the velocity, but this will not have a 
significant effect. The latter assumption is justified because 
an azimuthal velocity component in the post-shock gas will 
require one to be present in the pre-shock gas, and this will 
not be the case in a radially outflowing wind. 

With these assumptions, the centroid shift (v) and ve- 
locity range (f ma x — ti m i n ) of an emission line are given by 
l|Luhrj|l997l:IPollock et al]l2004D 

v — — vo cos (5 cos 7 (9) 

Smax - fmin = 2Av = 2v SU1 (3 SU1 7 (10) 

In particular, v comes from 'Term 1' in lLiihrsI ' equation (9), 
where his equals our ^ and his «wa equals our vo, while 
Av is equal to v* in his equ ation ( 8 ) [see also his equa- 
tions (10) and (11)]. Note that lLuhra Jl997f) does not use 7 
in his analysis, but instead expresses his velocities in terms 
of the orbital phase angle and inclination explicitly. 

If we identify Av in equation JTJ]) with the FWHM, 
then given the mean observed line shift and FWHM (y = 
-64 ± 12 km s" 1 , FWHM = 1240 ± 30 km s" 1 ; see Sec- 
tion H3J and 7 = 44° (see Section 03. we obtain j3 = 87° 
and wo = 1800 km s _1 . This inferred value of vq is rather 
large, given that the terminal velocity of the WR star wind 
is 1500 km s _1 and the O star wind is likely to be collid- 
ing well below its terminal velocity of 2300-2400 km s _1 . 
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Figure 12. Geometrical model of the wind-wind collision in 
7 2 Vel. The solid circles represent the two stars. The cone (with 
opening half-angle /3) represents the wind-wind interaction region 
(along which X-ray-emitting material is streaming at speed vq). 
The viewing angle 7 is the angle between the line of sight and the 
line of centres. 

However, the inferred value of vo is sensitive to the relation- 
ship between (i> ma x — Wmin) and the FWHM [we have used 
FWHM = 0.5(?; max — u m m) in a rather ad hoc manner to 
characterize the line widths], and hence vo may in fact be 
somewhat lower. This is not a problem for f3 because it is 
very insensitive to the relationship between (n ma x — Vmin) 
and the FWHM. 

Note that the mean observed line shift is compara- 
ble to the absolute wavelength accuracy of the HETGS 
(~100 km s _1 ; Chandra Proposers' Observatory Guide, Sec- 
tion 8.2.4). This means that the measured value of v (and 
hence the derived values of j3 and «o) may not be completely 
trustworthy. However, we can say that \v\ < 100 km s _1 
(since we would unambiguously be able to detect a larger 
value of v). This implies that /3 > 85°. The consequences of 
this are discussed later. 



5 LINE PROFILE MODELLING 

In an attempt to better understand the line shifts and widths 
observed in the Chandra spectrum of 7 2 Vel, we have calcu- 
lated model X-ray l ine profiles using the model described in 
iHenlev et all (I2003I) . This model uses data from a hydrody- 
namical simulation of the wind- wind collision. The simula- 
tions are two-dimensional and assume cylindrical symmetry 
about the line of centres. The stars' winds are assumed to be 
spherically symmetric, and orbital effects are assumed to be 
negligible (this latter assumption is discussed later). The line 
profiles are calculated by assuming each grid cell produces 
a thermally Doppler-broadened (i.e. Gaussian) line profile, 
the amplitude of which depends on the density and temper- 
ature of the gas in the cell. The centre of the Gaussian is 
shifted according to the line-of-sight velocity of gas. Since 
the simulations assume cylindrical symmetry, each grid cell 
actually represents a ring of gas. The line-of-sight velocity 
of the gas varies as a function of position around the ring. 
This is taken into account by dividing the ring into a number 
of sections (typically 200), and calculating the contributions 



from each section (the line-of-sight velocity of a given section 
depends on the 2-D velocity components in the appropriate 
hydro grid cell, the viewing angle 7, and the azimuthal angle 
(j) of the section around the ring). Absorption by the cold, 
unshocked winds of the stars is included in the calculations. 
For each emitting ring section, the integral of the density 
along the line of sight is calculated (where possible by tak- 
ing the density directly from the hydro grid; off the grid the 
density is extrapolated by assuming spherically symmetric 
winds). This is multiplied by the opacity (calculated for a 
solar abundance plasma using xstar 6 ) to give the optical 
depth r, and then the amplitude of the Gaussian profile pro- 
duced by that ring section is multiplied by e _T . The overall 
line profile is simply the sum of the contributions from the 
whole grid (summing over z, r and (f>). 

The major factors which affect X-ray line profiles 
from colliding wind binaries are the mass-loss rates, the 
wind speeds, the separation and the viewing angle (i.e. 
the angle between the line of sight and the line of cen- 
tres). As the orbit is well d etermined iSchmutz et aljFl997l : 
iDe Marco fc Schmutzll999l'l the separation (0.92-0.94 A.U.) 
and the viewing angle (7 = 44° ; see Section 14.51 can be 
specified a priori. The appropriate set of wind parame- 
ters to use in the hydro dynamical simulation is less cer- 
tain. iBarlow et alJ lll988T) measured the terminal velocity 
of the WR star to be 1520 km s _1 from th e half-width of 
the fo rbidden 12.8^m [Nell] emission line. ISt-Louis et alJ 
( 1993) inferred a very similar value (1550 km s _1 ) from their 
analysis of phase-varying absorption in IUE spectra. Given 
that the WR star wind is more powerful than the O star 
wind, the wind-wind collision will be close to the O star. 
As the separation is ~60i?» (WR) , the WR star wind is very 
likely to be at its terminal velocity at the wind-wind colli- 
sion. Conversely, the O star wind (vqq » 2300-2400 k m s" 1 ; 
IPrinia. Barlow fe Howarthll990l : lst-Louis et aHl993ft is un- 
likely to be at its terminal velocity at the wind-wind colli- 
sion, as evidenced by the disappearance of the high- velocity 
blue absorption wing from ultraviolet P Cygni profiles ob- 
tained by IUE when the O star is behi n d the WR star 
iSt-Louis e t al. 1993|) and bv lWillis et alt s 1119951) compar- 
ison of model X-ray spectra with their observed ROSAT 
spectra. Note that as the separation varies throughout the 

orbit, so too will the pre-shock velocities. 

As already mentioned in Section IStevens et alJ 

(1996) derived from ASCA data a mass-loss rate for the 
WR star of 3 x 10 ~ 5 M fl yr _1 , a factor of three lower 
than iBar low et all 's il988f) radio-det e rmine d value of 8.8 x 
10 -5 Mq yr _1 , which IStevens et al.l il996l) suggested may 
be due to the fact that radio observations tend to over- 
estimate mass-loss rates when the wind is clumped. It 
should be noted that both these analys es used larger, pre - 
Hipparcos d istances to j 2 Vel o f 460 pc iBarlow et al.lll988l) 
and 450 pc iStevens et al Il996ft . Since the radio-determined 
mass-loss rate scales a s D 3 ^ 2 , where D is the distance 
JWright fc BarlowHl975Tl. using the Hipparcos dis ta nce o f 
258 pc JSchaerer et alJll997h reduces IBarlow et alJ 's (Il988l) 
value to 3.7 x 10" 5 M Q yr" 1 . If one assumes the observed 
X-ray flux scales as M 2 , the mass-loss rate determined 
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using IStevens et alt s lll99rJ) method scales as D, and so 
their value reduces to 1.7 x 10 -5 Mq yr -1 . More recent val- 
ues for the mass-loss rate of the WR star (which use the 
Hipparcos distance) include 1.08 or 3.06 x 10 -5 Mq yr -1 
based on rad io data (assuming a clumped or smooth wind, 
respectively: iNueds. Crowther fe Willis] Il998l) . and 0.93 or 
3.3 x 10 -5 Mq yr -1 based on stellar atmosphere modelling 
of He I and Hell lines (the different valu es coming from 
differ ent codes and different filling factors; iDe Marco et alj 
2000). However, as noted in Section f3. II the Hipparcos dis- 
tance has recently been thrown into doubt by the discovery 
of an association of low-mass, pre-main sequence stars in 
the direction of 7 2 Vel, and the distance to Vel may be 
between 360 and 4 90 pc ijPozzo et aljEoOch . 

For the O star. lSt-Louis et alJ <ll993T) adopted a typical 
mass-loss rate for a 091 star (the t hen accepted spectral clas - 
sification) of 1. 3 x 10~ 6 M^ vr" 1 jHowaith fc Prini Jl98Sft . 
More recently. IDe Marco fc Schmut j Jl^gglT liave derived a 
much lower Mo of 1.8 x 10 -7 Mq yr -1 from hydrodynamical 
calculations of the radiatively driven wind. 

In hydrodynamical simulations of the wind-wind col- 
lision in 7 2 Vel, one must take into account radia- 
tive cooling. This is because the cooling parameter x 
(the ratio of the cooling timesc ale to the flow timescale; 
IStevens. Blondin fc Pollocklll992|) is less than unity for the 
WR star wind. One effect of this is that the wind- wind in- 
terac tion region can become highly unstable iStevens et alJ 
1992), leading to strong time- variability of the calculated 
line profiles. These instabilities may be partially numerical 
in origin, and so may lead to unphysical results. Further- 
more, if the cooling is very strong, which will be the case for 
a large ad opted mass-loss rate or low wind speed [see equa- 
tion (8) in lStevens et alJl992|| . the shocked gas will collapse 
into a thin sheet and will not be properly r e solved on the 
grid. This means the model of IHenlev et alJ 11200311 cannot 
be used, as the profiles are calculated after the hydrodynam- 
ical simulation and so do not take into account the energy 
that was radiated during the simulation. 

The line profile calculations presented here are based 
on simulations ca rried out with the hydrodynamical code 
cobr a (see, e.g., iFalle fc Komissarovl Il996t iPittard et alJ 
2003), which is second order accurate in space and time. Ra- 
diative cooling is implemented in the simulations. A small 
amount of numerical viscosity is used to damp the growth 
of a numerical instability which appears when shocks are 
stationary on the grid. We adopt mass-loss rates of Mwh 
1 x 10~ 5 Mq yr -1 and Mo = 5 x 10~ 7 Mq yr -1 , and wind 
speeds of 1500 km s _1 for both winds. For simplicity, we as- 
sume the winds are not accelerating, which is why we adopt a 
wind speed for the O star much lower than its terminal veloc- 
ity. Fig. ll3l shows density, temperature and speeds map from 
the resulting simulation. As can be seen, the cooling of the 
post-shock gas is resolved. Our simulations show little sign of 
instabilities or post-shock mixing and turbulence. The lack 
of any instability is unsurprising for two reasons. Firstly, the 
wind speeds are equal, so there will be very little slip along 
the contact discontinuity and the Kelvin-Helmholtz insta- 
bility will not arise. Secondly, the cooling time is long (espe- 
cially for the O star wind, which remains almost adiabatic), 
so the shocked gas does not collapse into a thin sheet and 
the thin-shell instability is not excited. Even if turbulence 
is present in reality, in order to affect the results discussed 
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Figure 13. Density (top), temperature (middle) and speed (bot- 
tom) maps from a hydrodynamical simulation of ■y 2 Vel. The WR 
star is in the bottom left-hand corner of the grid; the O star 
is half-way along the bottom axis. The adopted wind parame- 
ters are A/ WR = 1 X 10~ 5 Mq yr -1 , M = 5 X 10~ 7 Mq yr -1 , 
«WH = vq = 1500 km s . 



below the turbulent velocity would have to be a significant 
fraction of the bulk flow speed, which is unlikely to be the 
case. 

We have calculated profiles for the Lya lines and He- 
like resonance lines of Ne, Mg, Si and S. The line emissiv- 
ity is assumed to vary just with temperature; we do not, 
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for example, include the en hancement of the H e-like reso- 
nance line at high densities foorauet et al.ll200ll) . as it will 
not be important at the densities present in the simula- 
tion results. The line emissivities were taken from ATOMDB 
vl.1.0, which assumes solar abundances. To take into ac- 
count the different compositions of the winds, we normal- 
ize the emission from each wind using the product of the 
number density of the element in question and the num- 
ber density of electrons. The elemental number densities are 
calculated from the gas density p using the elemental mass 
fractions. For the Q star wind we assume solar abundances 
l)Anders fc Grev cssc 1989), while for the WR wind the mass 
fractions are calculated us i ng N(R) = 0, N(C){N(Ue ) 
= 0.14 jMorris et alJ ll99St ISchmutz fc De Marcc *" 
AT(N)/iV(He) = 1.0 x 10~ 4 folovd fc Sticklanc 
N(G)/N( 0) = 5 foe Marco et al.ll200CT) . MNel/MHe) = 
3.5 x 10~ 3 foessart et all2000h and N{S) /iV(He) = 6 x 10" 5 
foessart et alJl2000|) . [In the preceding. N(X) denotes the 
number abundance of element X.] The relative abundances 
of elem ents more massive than n eon are kept at their solar 
values l|Anders fc Grevesselll989T) . Note that we do not use 
elemental abundances measured in Sectional as the different 
spectral models used give significantly different abundances 
for some elements. However, in this case the shapes of the 
calculated profiles are not very sensitive to the abundances 
used. The electron density n c is calculated from the gas den- 
sity using p — 2mnn c /(l + X), where mn is the mass of 
a hydrogen atom and X is the mass fraction of hydrogen 
foowers fc Deemindll984l equations 2.22 and 2.23). Since 
the line fitting described in Section^Jtakes into account the 
presence of the fainter component of each Lyman line and 
the intercombination and forbidden lines of the fir triplets, 
the lines we calculate are singlets, enabling a direct compar- 
ison with the shifts and widths measured from the HETGS 
spectrum. 

Examples of line profiles calculated for 7 = 44° are 
shown in Fig. 1141 which shows the Si xm resonance line and 
the Si xiv Lya line. The lines are redward- skewed, rather 
than Gaussian. This skewing is intrinsic and due to the ge- 
ometry of the emitting region, rather than being due to ab- 
sorption, which does not have a strong effect on these lines. 
The lines are broad, each with a FWHM of ~1000 km s _1 , in 
good agreement with the widths measured from the HETGS 
spectrum. However, the lines are noticeably blueshifted, 
with centroid shifts of ~— 300 km s _1 , in contrast to the 
essentially unshifted lines in the observed spectrum. The 
calculated line blueshift decreases as 7 increases, but given 
the opening angle of the shocked region in Fig. 1131 increas- 
ing 7 much above ~40° means the X-ray emitting plasma 
would be observed through the WR star wind instead of 
the O star wind, and this is unlikely given the amount of 
absorption observed in the spectrum (the column density 
measured in Section [3] would probably be much larger if the 
wind-wind collision region were being observed through the 
WR star wind). 

Because the line profile calculations described here are 
based on 2-D hydrodynamical simulations, they do not take 
into account the effects of orbital motion on the shape 
of the wind-wind collision region. However, one can es- 
timate the importance of these effects by assuming the 
Coriolis force deflects the whole shock cone by some angle 
£ = tan -1 (i; or b/vwind) in the orbital plane, where v OI b and 



Wwlnd are the orbital and wind speeds, respectively. At the 
time of the Chandra observation, the speed of the O star in 
the rest frame of the WR star was ~200 km s _1 . If we as- 
sume Uwmd ~ 1000 km s" 1 , we find £ « 11°. This means \& 
[see equation ©] is decreased to 25°, and the viewing angle 
7 decreases to 36° . This actually increases the blueshift that 
the model predicts to ~350 km s" 1 (in this case the Coriolis 
force tends to deflect the shock cone such that there is more 
gas travelling towards the observer, increasing the predicted 
blueshift). Using 7 = 36° instead of 44° in Section [4.51 still 
predicts a very large shock opening half-angle (/3 > 86°). 
This simple estimate of the size of the effect of the Coriolis 
force does not take into account the fact that the wind-wind 
interaction region will change shape, rather than simply be- 
ing deflected en masse (the whol e region becomes be nt round 
into an Archimedean spiral; e.g. iTuthill et al]ll999l) . Never- 
theless, it appears that Coriolis effects cannot explain the 
discrepancy between the observed and calculated line shifts. 

The discrepancy between the observed and calculated 
line shifts is most likely due to the opening half-angle of the 
shocked region being ~40° in our hydrodynamical simula- 
tion (see Fig. 1131 , whereas the simple geometrical model of 
the emission lines in Section [4.51 implies a shock-cone open- 
ing half-angle of > 85°. In a CWB with non-accelerating 
winds, /3 is a function only of the wind momentum ratio 
(MwrVwr/A/o^o). The shock opening angle in our hy- 
drodynamical simulation therefore depends on the adopted 
values of the wind parameters. However, the shock open- 
ing half-angle inferred from our simple geometrical model 
(P > 85°) implies approximately equal wind momenta, 
which is not the case for any sensible set of (constant veloc- 
ity) wind parameters. We therefore suggest that this large 
implied openi ng angle may be evidence of sud den radia- 
tive braking iGavlev. Owocki fc CranmeJ IT997). in which 
the wind of the WR star is rapidly deceler ated when it 
encou nters the radiation field of the O star. iGavlev et alJ 
lll997l) suggest that the wind-wind collision in -y 2 Vel could 
be significantly affected by sudden radiative braking, as they 
find that the wider separation than in V444 Cygni (the main 
system they study) offers greater opportunities for braking. 



6 DISCUSSION 

The broad-band X-ray spectrum of Vel indicates that a 
wide range of temperatures is present in the X-ray-emitting 
plasma, covering a range of ~4-40 MK. The temperature 
distribution below ~4 MK is poorly contrained, because of 
the low number of counts above ~14 A. The differential 
emission measure (Fig. indicates the temperature distri- 
bution falls off above ~40 MK. This temperature range im- 
plies shock speeds of 500-1700 km s" 1 assuming solar abun- 
dances, or 300-1100 km s _1 assuming WC8 abundances. 
The upper values of these ranges are reasonable given the 
discussion of wind speeds in the previous section. The post- 
shock temperature corresponding to Uoo(O) = 2400 km s _1 
is 80 MK. The lack of any evidence for gas at such high tem- 
peratures suggests that the O star wind is not colliding at 
its terminal velocity near the line of centres [the shocked gas 
near the line of centres is the densest in the collision region, 
and so if there is any very hot (~80 MK) gas it should make 
a significant contribution to the X-ray emission]. However, 
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Figure 14. Model Si XIV Lya (left) and Sixm resonance (right) 

it should be noted that the sensitivity of the HETGS falls 
off at short wavelengths (A < 2 A) and so it is less sensitive 
to emission from very hot gas. To investigate this further, 
we synthesized a spectrum by adding a 80 MK apec compo- 
nent to our best-fitting 2T apec model (model B in Tablc0. 
With its emission measure [expressed in terms of the helium 
density (EM = J n c nn c dV), because of the lack of hydro- 
gen in the WR star wind] fixed at 2 x 10 53 cm -3 , this model 
leads to significant Fexxv emission near 1.85 A which is not 
observed. This implies that the emission measure of the very 
hot (80 MK) gas is at least a factor of 5 smaller than those 
of the two apec components in Tabled 

The lower limits of the implied shock speeds do not 
imply low wind speeds per se, but instead that parts of 
the wind are being shocked obliquely as is expected from 
the geometry of the wind-wind collision region. This oblique 
shocking also explains (at least qualitatively) the shape of 
the DEM in Fig. |K| If all the gas were being shock-heated 
to ~40 MK at the shock apex and then allowed to cool ra- 
diatively, one would expect the DEM to have a shape that 
is the inverse of the cooling function, i.e., at temperatures 
where the cooling function is low, there will be more gas 
(as it will take longer for gas to cool through this temper- 
ature) and so a larger emission measure (and vice versa). 
Therefore, one would expect the emission measure to de- 
cr ease from ^20 MK to lower tempe ratures (see, e.g., fig. 3 
in lAntokhin. Owocki fc Brownll2004) . Instead, the observed 
emission measure rises from ~20 MK to ~8 MK. This means 
that gas is being fed in at lower temperatures, which is con- 
sistent with the oblique shocking that occurs away from the 
line of centres. 

The column density derived from the broad-band fitting 
(N h ~ 2 x 10 2 2 cm" 2 ) is consistent with that measured 
by iRauw et alJ i2000T) from an ASCA observation at the 
same phase (N K = 2.4 x 10 22 cm" 2 at $ = 0.078). The 
measured value indicates the wind-wind interaction region 
is being observed through the wind of the O star, which is 
what is expected given the viewing orientation. 

As has already been stated, the broad-band spectral 
models used in Sectional give poor fits to the data (xt ~ 2). 
Furthermore, they can at best only characterize the general 
properties of the X-ray-emitting plasma. A better approach 
is to calculate the X-ray emission expected from the hydro- 
dynamics of the wind-wind collision. This may be done using 
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nume rical hydrodynamical simulations (e.g. IStevens et al.l 
1996), in a similar fashion to the method used in Section |S] 
An alternative is to use th e semi-analytical approach devised 
bv lAntokhin et all |2004) . Firstly, the shape of the interac- 
tion region is calculated using wind momentum balance. The 
shocks are assumed to be highly radiative, so the shocked 
material collapses into a thin layer and the two shocks and 
the contact discontinuity (CD) are treated as effectively syn- 
onymous. At the shocks (which are treated as locally pla- 
nar), the wind is shock-heated to some temperature T s and 
then cools rapidly (effectively to T — 0) in a thin layer 
between the shock and the CD; i.e., all the kinetic energy 
associated with the velocity component locally perpendic- 
ular to the shock is radiated away. The former approach 
works best for adiabatic shocks (the instabilities of radia- 
tive shocks lead to substantial mixing between hot and cold 
material, the physically appropriat e degree of whic h is dif - 
ficult to predict a priori), whereas lAntokhin et alt s (l2004h 
method assumes fully radiative shocks. In 7 2 Vel, the O star 
wind is effectively adiabatic, and while the WR star wind is 
radiative, it is not highly radiative (the post-shock cooling 
is resolved in Fig. 1131 . Nevertheless, a useful next step in 
our understanding will be the application of both methods 
to the Chandra HETGS spectrum of 7 2 Vel. Since veloc- 
ity information is available when synthesizing spectra using 
either method, it should also be possible to calculate line 
shapes self-consistently. Both methods may be used to gen- 
erate grids of spectra for a range of parameters (e.g. M, Woo). 
These can then be fit to the spectra in order to constrain 
wind parameters. In order to fit to the wealth of detail in the 
HETGS spectrum, the most up-to-date atomic data avail- 
able must be used to synthesize the spectra. However, such 
an analysis is beyond the scope of the current paper. 

From the emission lines in the X-ray spectrum, we mea- 
sured temperatures by comparing the fluxes of lines from H- 
like and He- like ions of a given element, and by using G ratios 
derived from He- like fir triplets. For Ne and Mg, the temper- 
atures measured by the two methods are in good agreement 
with each other. However, the Sixiv:Sixill flux ratio implies 
a significantly higher temperature (12 MK) than the Sixm 
G ratio (5 MK). As already stated, this may simply indi- 
cate that the Si XHI [ionization potential (LP.) = 2.438 keV; 
lDappenl2 000l originates from somewhat cooler regions than 
the Si xiv emission (LP. = 2.673 keV). However, it could 
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also be evidence of non-equilibrium ionization (NEI). This 
could tie in with the evidence that the Mgxi emission (LP. 
= 1.7618 keV) seems to originate closer to the O star than 
the Si xiii emission, as evidenced by the R ratios derived 
from their He-like fir triplets. Furthermore, the G ratios in- 
dicate that the Mgxi emission originates from hotter gas 
(T e = 7 MK) than the Si XIII emission (T e = 5 MK). This 
implies it originates nearer the line of centres (and hence 
nearer the O star), as the temperature decreases monoton- 
ically away from the line of centres (see Fig. 1131 . This is 
counter to what is expected if collisional ionization equilib- 
rium holds, as in that case the Mgxi emission would origi- 
nate in cooler gas further from the line of centres than the 
Si xm emission. 

One can assess whether or not NEI is important in a 
colliding wind binary by deriving an ionizat ion parameter 
V ana logous to the cooling parameter x of IStevens et alJ 
(1992). In the adiabatic limit, NEI effects will be impor- 
tant if the ionization equilibration timescale is larger than 
the flow timescale, whereas in the radiative limit they will 
be important if the ionization equilibration timescale is 
larger than the cooling timesc ale. The ioni zation equilibra- 
tion timescale t\ CCL is given by jMasailll994ft 



^Ieq — 
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where n c is the electron number density. For the flow 
timescale ta ov , we shall use 
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where d is the distance from the star to the c ontact discon- 
tinuit y and c s is the post-shock sound speed iStevens et alJ 
|l992j). The cooling t imescale t C ooi is given approximately by 
JStevens et"aflll992ft 
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where T is the post-shock temperature, n is the number 
density of the wind at the shock and A is the cooling function 
(~ 2 x 1 0~ 23 erg s -1 cm 3 for solar-abundance material; see 
fig. 10 in lStevens et al]ll992ft . 

Note that ifl ow in equation 112H is the timescale for the 
post-shock gas to travel a distance ~d away from the line 
of centres, whereas the temperature in the post-shock gas 
may decrease away from the line of centres such that the 
ionization balance for a given element is expected to go from 
being dominated by H-like ions to being dominated by He- 
like ions in a distance less than d. It is therefore possible 
that the above flow timescale is too long to be appropriate. 
Nevertheless, with this choice of tg ow one finds that in the 
adiabatic limit 
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whereas in the radiative limit 
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where fi c is the average mass per electron in a.m.u. (1.2 
for solar abundance material and 2 for WC material), /i is 
average mass per particle in a.m.u. (0.61 for solar abundance 
material and 1.5 for WC material), di2 is the distance from 



the star to the contact discontinuity in units of 10 cm, 
v$ is the wind speed in units of 1000 km s _1 , and M_7 is 
the mass-loss rate in units of 10 -7 Mq yr _1 . If ib <C 1 then 
collisional ionization equilibrium holds. 

To determine ip we assume Mwr = 3x 10 -5 Mq yr~\ 
Mo = 5 x 10 -7 Mq yr" 1 and vwr = «o = 1500 km s _1 , 
and calculate di2(WR) and di2(0) by assu ming ram pres- 
sure balance at the contact discontinuity iStevens et alJ 
Il992ft . Using equation (TEH , we find Vwr = 0.004 and 
ipo — 0.02, whereas equation I15H gives Vwr = 0.01 and 
ipo — 0.04. Note that the value of Vwr in the radiative limit 
is probably underestimated by a factor of a few, because 
A is larger for WC material than so lar abundance mate- 
rial (see fig. 10 in lStevens et al.lll992ft . The winds in Vel 
are neither adiabatic nor completely radiative, and so the 
true values of Vwr and Vo lie between the above values. 
This therefore suggests that collisional ionization equilib- 
rium should hold for both stars' winds in 7 2 Vel, in contrast 
to what we infer from the Mgxi and Sixm fir triplets. A 
detailed self-consistent calculation of the ionization balance 
in the wind- wind collision is needed to resolve this issue. 

As well as inferring the temperatures and locations of 
different regions of X-ray emitting plasma from the emission 
line spectrum, we have used the line shifts and widths to in- 
fer the geometry of the wind-wind collision. The essentially 
unshifted lines imply a large shock-cone opening half-angle 
f3 > 85° . As stated pre viously, this may be evidence of sud- 
den radiative braking iGavlev et al 1 11997ft . Note that the 
calculation that leads to this opening angle does not take 
into account absorption in the cool, unshocked wind of the 
O star, which can have a profound affect o n the shape of 
the X-ray emission lines (iHenlev et alj|2003ft . However, the 
calculations in Section |S] while they do not reproduce the 
observed line shifts, do indicate that absorption does not 
have a strong effect on the line profiles in j 2 Vel. 

Earlier observations of 7 2 Vel, however, do not im- 
ply such a wide shock opening half-angle. From the du- 
ration of t he period of enhan ced emission in the ROSAT 
lightcurve, IWillis et alJ |l995) infer a shock opening half- 
angle of ~25° . Their estimate implicitly assumes an inclina- 
tion of 90°, but even using the inclination we have adopted 
throughout this paper (i = 63°: lDe Marco fc Schmutzil999ft 
only increases the inferred half-angle to ~50°. Further- 
more, if the shock opening half-angle is as large as 85°, 
iRauw et all s {2000) ASCA observation at phase $ = 0.978 
would have seen the wind-wind collision through the O star 
wind, and their measured column density would have been 
similar to that which they measured for their <3? = 0.078 ob- 
servation (Nh — 2.4 x 10 22 cm -2 ), rather than three times 
larger (JV H = 7.5 x 10 22 cm -2 ). 

This discrepancy between the shock opening angle in- 
ferred from the line shifts and the previous X-ray obser- 
vations of 7 2 Vel leads one to speculate that the intrinsic 
emission of the O star is contributing a significant amount 
to the observed X-ray spectrum (the WR star is unlikely to 
contribute any intrinsic emission, as no sin gle WC star has 
ever been convincingly detected in X-ravs: IOskinova et alJ 
2003). The standard wind-shock model, with X-rays being 
emitted from shocks distributed throughout the wind due 
to instabilities in the line-driving force (e.g. lOwocki et alJ 
Il988ft. predicts broad, blueshifted, blu e-skewed line profiles 
Jlgnacel lioOll : lOwocki fc Cohenl l200lt) . as the red-shifted 
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emission from the far side of the wind is attenuated by 
the wind. Such profiles h ave been observed fro m the well- 
studied 04f star £ Pup (|Cassinelli et alJl200l|) . and have 
been successfully fitted with a simple model of wind absorp- 
tion llKramer. Cohen fc Owockil2003h , but in general a clear 
picture of X-ray line emission from early-type stars has yet 
to emerge. The models for X-ray line profiles from single O 
stars discussed above cannot be applied to y 2 Vel because 
they rely upon spherical symmetry, which is not the case 
in j 2 Vel because of the wind-wind collision. However, if 
the lines were coming from the O star wind, the FWHM 
(1200 km s _1 ) implies that most of the emission originates 
from where the line-of-sight velocity is between —600 and 
+600 km s _1 (i.e. from where \v\ <«oo/4). If one assumes 
a f3 velo city law of the form v(r) — ?;oo(l — R*/r)^ with 
= 0.8 JSt-Louis et al j ll99^ . this implies the line emission 
mostly originates less than 1.27?, (O) from the centre of the 
O star, which is inconsistent with the minimum radii of for- 
mation inferred from the 7? ratios of He-like ions. Of course, 
this argument says nothing about the lines from H-like ions, 
but given that the widths and shifts are similar for all ions 
(Fig. |HJ , it is reasonable to assume a common origin for all 
the lines. 

The I/x/iboi ratio could be used to assess the impor- 
tance of the O star's i ntrinsic emissio n to the total emission 
observed from j 2 Vel. iPollockl il987f) measured Lx/Lboi = 
0.44 x 10~ 7 for 7 2 Vel with Einstein, which is typical for lone 
O sta rs (Lx/Lboi ~ 10~ 7 : IChlebowskilll989l: iMoffat et ail 
2002). Model B in Table Q gives an intrinsic X-ray luminos- 
ity in the Einstein (0.2-3.5 ke V) band of 9.6 x 10 32 erg s" 1 . 
Using Lboi(O) = 2.1 x 10 5 L Q foe Marco fc Schmutdll999Tj 
gives Lx/Xbol = 12 x 10 -7 , which is much higher than the 
Einstein value and which seems to indicate that colliding 
wind emission dominates. Unfortunately, due to the low 
number of counts at low energies, the luminosity at lower 
energies is very poorly constrained, and so this conclusion 
may be unreliable. However, the high temperatures inferred 
from the spectrum (up to ~40 MK) are evidence of colliding 
wind emission, as the standard wind-shock model does not 
produce strong enough shocks to produce such high temper- 
atures. 



7 SUMMARY 

We have carried out a new analysis of an archived Chan- 
dra HETGS spectr um of 7 2 Vel. As w ith the previously 
published analysis JSkinner et alJ l200ll) we find the lines 
are essentially unshifted, with a mean FWHM of 1240 ± 
40 km s -1 . The lack of strong absorption means the line-of- 
sight is through the O star wind. We find a wide range of 
temperatures is present in the X-ray-emitting plasma, from 
~4 to ~40 MK, which is in good agreement with the range 
of temperatures expected from a hydrodynamical simula- 
tion of the wind-wind collision. However, line profile cal- 
culations based on such simulations imply the lines should 
be blueshifted by a few hundred km s _1 . We suggest that 
the unshifted lines may imply a wider shock opening half- 
angle (> 85°) than that which is seen in the hydrodynamical 
simulations. This may be evidence of sudden radiative brak- 
ing. However, such a wide shock opening half-angle seems 
to disagree with earlier ROSAT and ASCA observations of 



7 2 Vel. A full radiation hydrodynamics simulation of y 2 Vel 
is needed to investigate what effect radiative braking has 
on the morphology of the wind-wind collision and how the 
degree of this effect may change with orbital separation (see 
|Pittardlll99l). 

Unlike ISkinner et al.l (1200 lb . we find evidence that the 
Mgxi emission emerges from hotter gas closer to the O star 
than the Sixill emission, which may be evidence of non- 
equilibrium ionization. Future work will include trying to 
model the ionization balance in colliding wind binaries self- 
consistently. This may be important fo r the study of othe r 
collidi ng wind systems, s uch as WR140 dPollock et al.l2o'o4T) . 

iGavlev et alJ il997Tl point out that the eccentricity of 
7 2 Vel allows the degree of braking to be studied as a func- 
tion of the orbital separation. Unfortunately, the strong ab- 
sorption in the WR star wind means that a long (~200 ks) 
Chandra pointing would be needed to obtain a high-quality 
spectrum at phases other than when the O star is in front. 
However, 7 2 Vel is an ideal target for Astro-E2, due to its 
greater effective area at higher energies. We have proposed a 
75 ks Astro-E2 observation of 7 2 Vel when the O star is not 
in front (PI: D. Henley). This will enable us to measure line 
shifts and widths at a different orbital phase from the Chan- 
dra observation. From this we can investigate how the shock 
opening angle varies with orbital phase. By comparing this 
with the shock opening angle inferred from radiation hydro- 
dynamics simulations, we will be able to place constraints 
on the coupling between the O star radiation field and the 
WR star wind, which will in turn provide insights into the 
poorly understood physics of WR wind driving. 
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